#include "llcomp_llc.h" 

    

     
    
     void loopParallelFor0_support(  int *_i,  int *_j,  double *x,  double *_sum,  int *_n,  double **a,  double *b,  double *_resi,  double *new_x) ;
     void loopParallelFor0_support(  int *_i,  int *_j,  double *x,  double *_sum,  int *_n,  double **a,  double *b,  double *_resi,  double *new_x)  {
                 int i = *_i;
                 int j = *_j;
                 double sum = *_sum;
                 int n = *_n;
                 double resi = *_resi;
        /* Begin template: forall_init_unit.tmpl */

/* $0 => Variable del bucle
 * loop_initial => Valor inicial de la variable del bucle B
 * loop_finalice => Valor final de la variable del bucle C
 */

{
    int llc_meanT, llc_rnT;
    register int llc_i;

    llc_nT = ((n - 1) - (0) + 1);
    /* Medium number of tasks */
    llc_meanT = llc_nT / LLC_NUMPROCESSORS;
    /* Number of unassigned tasks */
    llc_rnT = llc_nT % LLC_NUMPROCESSORS;

    LLC_CHECK_MALLOC (llc_F, (LLC_NUMPROCESSORS + 1), int);

    llc_F[0] = 0;
    for (llc_i = 1; llc_i < (LLC_NUMPROCESSORS + 1); llc_i++) {
        llc_F[llc_i] = llc_F[llc_i-1] + llc_meanT + (((llc_i - 1) < llc_rnT) ? 1 : 0);
    }
}

/* End template: forall_init_unit.tmpl */
/* Begin template: forall_exec_unit.tmpl */

/* loop_var => Variable del bucle
 * loop_initial => Valor inicial de la variable del bucle
 * loop_stmt => Bloque del bucle
 */

{
    int llc_nP;
    int llc_grp_save;
    register int llc_i;

    llc_nP = LLC_NUMPROCESSORS;
    llc_grp_save = LLC_NAME;
    LLC_NUMPROCESSORS = 1;
    LLC_NAME = 0;
    
    i = (0) + llc_F[llc_grp_save];
    for (llc_i = 0; llc_i < LLC_PROCSGRP(llc_grp_save); llc_i++, i++) {
        
{

sum= 0.0;
  for (j= 0; (j < n); (j++ )) sum+= (a[i][j] * x[j]) ;
resi+= fabs((sum - b[i]) ) ;
sum+= ((-a[i][i] ) * x[i]);
new_x[i]= (((b[i] - sum)) / a[i][i]);
}

;
    }

    LLC_NAME = llc_grp_save;
    LLC_NUMPROCESSORS = llc_nP;
}

/* End template: forall_exec_unit.tmpl */
/* Begin template: forall_comm_unit.tmpl */

/* loop_var => Variable del bucle
 * loop_initial => Valor inicial de la variable del bucle
 * loop_finalice => Valor final de la variable del bucle
 * forall_comm_res_init => template dinamico que se utiliza para reservar memoria para las comunicaciones
 * forall_comm_res_pack => template dinamico que se utiliza para realizar las comunicaciones
 */

{
    register int llc_i, llc_j;
    long llc_buf_size = 0, llc_size_tmp = 0;
    char *llc_buf = NULL, *llc_buf_ptr = NULL;

    for (llc_i = 0; llc_i < LLC_NUMPROCESSORS; llc_i++) {
        llc_size_tmp = 0;
        for (llc_j = 0, i = (0) + llc_F[llc_i]; llc_j < LLC_PROCSGRP(llc_i); llc_j++, i++) {
            llc_size_tmp += (1) * sizeof ((&new_x[i])[0]);

        }
        llc_buf_size = LLC_MAX(llc_buf_size, llc_size_tmp);
    }
    
    if (llc_buf_size > 0) {
        
        LLC_CHECK_MALLOC (llc_buf, llc_buf_size, char);
        
        for (llc_i = 0; llc_i < LLC_NUMPROCESSORS; llc_i++) {
            if (llc_i == LLC_NAME) {
                    
                /* Simple element. Don't pack */
                if (LLC_PROCSGRP(llc_i) == 1) {
                    i = (0) + llc_F[llc_i];
                    MPI_Bcast(&new_x[i], ((1) * sizeof((&new_x[i])[0])), MPI_BYTE, llc_i, *llc_CurrentGroup); 
                }
                /* Multiple elements. Pack */
                else {
                    llc_buf_ptr = llc_buf;
                    
                    for (llc_j = 0, i = (0) + llc_F[llc_i]; llc_j < LLC_PROCSGRP(llc_i); llc_j++, i++) {
                        /* Packing &new_x[i]... */
memcpy (llc_buf_ptr, &new_x[i], (1) * sizeof ((&new_x[i])[0]));
llc_buf_ptr += (1) * sizeof ((&new_x[i])[0]);

                    }
                    
                    MPI_Bcast(llc_buf, (llc_buf_ptr - llc_buf), MPI_BYTE, llc_i, *llc_CurrentGroup);
                }
            }
            else {
                /* Simple element. Don't unpack */
                if (LLC_PROCSGRP(llc_i) == 1) {
                    i = (0) + llc_F[llc_i];
                    MPI_Bcast(&new_x[i], ((1) * sizeof((&new_x[i])[0])), MPI_BYTE, llc_i, *llc_CurrentGroup);
                }
                /* Multiple elements. Unpack */
                else {
                    llc_size_tmp = 0;
                    
                    for (llc_j = 0, i = (0) + llc_F[llc_i]; llc_j < LLC_PROCSGRP(llc_i); llc_j++, i++) {
                        llc_size_tmp += (1) * sizeof ((&new_x[i])[0]);

                    }
                    
                    MPI_Bcast(llc_buf, llc_size_tmp, MPI_BYTE, llc_i, *llc_CurrentGroup);
                    
                    llc_buf_ptr = llc_buf;
                    
                    for (llc_j = 0, i = (0) + llc_F[llc_i]; llc_j < LLC_PROCSGRP(llc_i); llc_j++, i++) {
                            /* Unpacking &new_x[i]... */
memcpy (&new_x[i], llc_buf_ptr, (1) * sizeof ((&new_x[i])[0]));
llc_buf_ptr += (1) * sizeof ((&new_x[i])[0]);


                    }
                }
            }
        }
        free (llc_buf);
    }
}

/* End template: forall_comm_unit.tmpl */
/* Begin template: forall_comm_red.tmpl */

{
    register int llc_i;
    int llc_buf_size = 0;
    char *llc_buf = NULL, *llc_buf_ptr = NULL, *llc_buf_tmp = NULL;
    MPI_Status llc_status;

    llc_buf_size += sizeof (double );
 

    if (llc_buf_size > 0) {
        
        LLC_CHECK_MALLOC (llc_buf, llc_buf_size, char);

        if (LLC_NAME == 0) {
            for (llc_i = 1; llc_i < LLC_NUMPROCESSORS; llc_i++) {
                MPI_Recv (llc_buf, llc_buf_size, MPI_BYTE, llc_i, LLC_TAG_REDUCE_DATA, *llc_CurrentGroup, &llc_status);
                llc_buf_ptr = llc_buf;

resi += (*(double  *) llc_buf_ptr);
llc_buf_ptr +=  sizeof(double );
 
            }

            llc_buf_ptr = llc_buf;

memcpy (llc_buf_ptr, &resi, sizeof(double ));
llc_buf_ptr += sizeof(double );
 
            
            MPI_Bcast (llc_buf, llc_buf_size, MPI_BYTE, 0, *llc_CurrentGroup);
        }
        else {
            llc_buf_ptr = llc_buf;

memcpy (llc_buf_ptr, &resi, sizeof(double ));
llc_buf_ptr += sizeof(double );
 
            
            MPI_Send (llc_buf, llc_buf_size, MPI_BYTE, 0, LLC_TAG_REDUCE_DATA, *llc_CurrentGroup);
            MPI_Bcast (llc_buf, llc_buf_size, MPI_BYTE, 0, *llc_CurrentGroup);
            
            llc_buf_ptr = llc_buf;

memcpy (&resi, llc_buf_ptr, sizeof(double ));
llc_buf_ptr += sizeof(double );

        }	

        free (llc_buf);
        free (llc_buf_tmp);
    }
}

/* End template: forall_comm_red.tmpl */
/* Begin template: forall_exit.tmpl */

    free(llc_F);
    llc_F = NULL;

/* End template: forall_exit.tmpl */

                *_i = i;
                *_j = j;
                *_sum = sum;
                *_n = n;
                *_resi = resi;

    }

    